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Abstract 

For various two dimensional non linear a models, we present a direct com- 
parison between the (3 functions computed with the 2 + e renormalization 
group and the (3 functions measured by Monte Carlo simulations. The the- 
oretical and measured (3 functions match each other nicely for models with 
a trivial topology, yet they disagree clearly for models containing topological 
defects. In these later cases, they are compatible with a phase transition at 
a finite temperature. This indicates that the global properties of the mani- 
fold do matter, in contradiction with the assumption used in the 2 + e RG 
computation. 
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In this letter, we want to address problems related to phase transitions in two dimensions 
(d = 2) (for models with continuous symmetries), the various non linear cr models (NLcrM), 
the renormalization group (RG), and the possible role of topological defects. First, we will 
review carefully our theoretical understanding of the various points just mentioned. Then, 
we will present extensive Monte Carlo simulations for different NLcrMs, and in particular a 
direct comparison with the RG j3 function at d — 2. 

The most studied models of phase transitions with continuous symmetry are the S^ 1 
Heisenberg models (or the n-components (p 4 models, as they are widely believed to be in 
the same universality class). In two dimensions, the accepted picture is that they have a 
phase transition at T c = 0, with a correlation length diverging as £ ~ exp(c/T) (we will 
discuss below the special case of the S 1 or X-Y model). This picture has been extended to 
essentially every two dimensional models with continuous symmetry. Of particular interest 
are the non linear cr models. In a NLcrM, the field has values in a manifold M = G/H, taken 
as the coset of a Lie group G by a subgroup H. These models are constructed so as to have 
a phase transition from a high temperature phase with symmetry G broken down to a low 
temperature phase with symmetry H. Their importance comes from the fact that, if there 
is universality in this larger class of models, the NLcrM is the canonical model to describe a 
phase transition with symmetry G broken down to H. The NLcrMs correspond to the limit 
where all massive modes are neglected. The simplest example of a NLcrM is constructed 
with G = 0{n) broken down to H = 0(n - 1), the manifold M = 0(n)/0(n - 1) = S^ 1 
is the n — 1 dimensional sphere. A point on a unit sphere is conveniently described by a 
n-components unit vector, and we find again the usual S 71-1 Heisenberg model. 

Polyakov initiated the RG computations in d = 2 + e for the Heisenberg models. 
The method was then further applied to all NLcrMs 0. The basic idea is to expand the 
Hamiltonian in a spin waves perturbation in the low temperature phase, i.e. to expand in 
the small displacements around an ordered state. Clearly, this method probes only the local 
structure of the manifolds. This is embodied by the resulting /3 function, which depends 
only on the curvature of the manifolds (at two loops order), a purely local geometrical 
quantity. In two dimensions, the computations predict generically a phase transition at 
T c = 0, with a correlation length diverging exponentially fast. However, the basic idea may 
appear dubious, because the transition between a low temperature phase (which depends 
only on the local properties of M) and a high temperature phase (which depends on the 
global properties of M) ought not to depend only on the local properties of the manifold. 
Strictly in two dimensions, the low temperature phase is even nonexistent, and yet, using 
this fictitious phase, we obtain predictions in the high temperature phase! The importance 
of the 2 + e RG computations comes from the fact that, except for the Heisenberg case, we 
do not have that many other reliable methods to obtain information about phase transitions 
(see e.g. |§ for a discussion of the problems related to the computations from the upper 
critical dimension and for the large n-expansion for the Stiefel model). To make it worse, 
we are often interested in three dimensional systems, or in quantum systems, i.e. to take 
e = 1 in the d = 2 + e expansion. Therefore, we need to have a reliable understanding of 
what happens at e = 0. Finally, let us note that, if the perturbative RG gives a prediction 
for the asymptotic regime (i.e. £ — » oo), it gives no estimate for the correlation length (or 
temperature) at which the asymptotic regime is reached. 

The role of topological defects in two dimensions is a related lingering problem. Most 
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of our comprehension comes from the X-Y or S 1 Heisenberg model. This model is unique 
for several reasons. The main one is that the manifold is of dimension one, and therefore 
possesses no curvature. This absence of curvature arises because the group G ~ 0(2) is 
Abelian. The related spin waves (3 function is then identically zero. Yet, this model possesses 
defects as Ili(5' 1 ) = Z2, and this drives the Berezinskii-Kosterlitz-Thouless transition [|J. 
Another peculiarity of the X-Y model is the Villain approximation ||, which decouples 
spin waves and defects, allowing to build a picture of the phase transition in terms of the 
defects only. Other NLcrMs have topological defects, as the related manifolds have some 
non trivial homotopy groups. Yet all models other than the X-Y model have curvatures. It 
has not been possible, up to now, to build an effective Hamiltonian for the spin waves and 
the topological defects (with their interactions) starting from the bare models. However, 
according to the RG argument, the spin waves alone are enough to disorder the system, and 
the topological defects are simply irrelevant. Thus, the NLuMs offer a larger testing ground 
for two dimensional phase transitions, the RG, and for the role of topological defects, away 
from all the peculiarities of the S 1 Heisenberg model. 

Essentially, the only tool at our disposal to test the various theoretical ideas presented 
above is the Monte Carlo (MC) simulation. In this letter, we present a direct comparison 
between the (3 function computed with the RG and the results of MC simulations. The 
basic idea is quite simple. Let us consider a model with cut-off A and one coupling constant, 
called T. The correlation length £(T, A) obeys the RG equation (A<9a + /3t<9t)£ = 0. As the 
correlation length is of dimension -1 we have A<9a£ = — £ and the previous relation can be 
put into the form (3t = l/$rln£. Then, the correlation length can be measured in a MC 
simulation, from which we can compute l/#rln£ and compare with the RG (3t function. 
Let us emphasize that there are no adjustable parameters in this process. The difficulty 
lays in the MC simulations because the correlation length is a poorly self averaging quantity, 
from which we have to take the logarithm and then a derivative. This will amplify the noise, 
and therefore we will need to have outstanding statistics in the MC simulations. 

Having presented the basic idea, let us turn to the general case. For a model defined with 
a cut-off A and a set of coupling constants g = the (bare) n points vertex functions 
satisfy a Callan-Zymanzik equation 



n 

A<9 A + (3id gi - -77 



rW(k j]9 ,A)~o (i) 



with (3i(g) = Ad\gi, T](g) = Ad\\nZ, and Z is used to renormalize the field. The right hand 
side is not strictly zero but it depends on a power of the cut-off, and is neglected according 
to a text book argument. This relation can be put into the form -77 = [A<9a + Pid gi ) lnT^. 
From the bare two points vertex function, we can compute the susceptibility x' 1 — T( 2 \k = 
0) and the correlation length £ 2 = d^p. In r^ 2 ^ . Combining these equations, we obtain 
a (Callan-Zymanzik) equation for the correlation length [A<9a + A<9 9i ]£ = and for the 
susceptibility [Ad A + ftid gi + f]]x — 0. The operator A<9a measures the dimension of the 
object it is applied on, and Ad\\n^ = —1, A^lnx = —2. Inserting in the above equation, 
we obtain our starting relation 

Wn£ = l (2) 
(3id g . lnx = 2 - r){g) 



3 



For a model with one coupling constant T, this becomes 



(3 T = l/d T \Tii (3) 



For a model with two coupling constants T and a, we obtain for the correlation length 

ffr _ 1 
1 — /3 a ^a In £ 9t In £ 



(5) 



In these relations, the left hand side can be computed with the RG, and the right hand side 
can be measured by MC simulations. 

Let us now define the NLcrMs we are working with. The space is a square lattice of size 
L 2 , with periodic boundary conditions. On each site x of the lattice, there is a field <p(x) 
with values in a manifold. The first set of models is the Stiefel model V ntP , corresponding to 
the group G = O(n), the subgroup H = 0(n — p) and the manifold V n>p — 0(n)/0(n — p) 
||. The field if is a n x p matrix, with the constraint ip + <f = Hp. The field can also be 
viewed as p orthonormal n-components vectors or as a p frame in IR™. The Hamiltonian 
is 

H = - £ Tr^ + (x) <p(y) = - £ X^iO) - <Pi(v) ( 6 ) 
<x,y> <%,y> i 

and the sum is over each link of the lattice. For p — 1, this is just the n-components 
Heisenberg model 5 n_1 ~ V n>1 . Let us note the special cases V^ n _i ~ SO(n) and V n>n ~ 
0(n). We will also use the SO (3) chiral model, corresponding to n = p = 3, but where 
the field obeys the supplementary constraint detip = +1. The second set of models is 
the Grassmann HP™ -1 model, corresponding to the group G = 0(n) and the subgroup 
H = 0(1) x 0(n — 1). The field is a n-components unit vector, and the Hamiltonian is 

h = - £ (*>(*) • v(y)f (7) 

<x,y> 

We made extensive MC simulations for all models with a manifold of dimension D — 1, 
namely V^i ~ 5* 1 (or X-Y model), ~ 0(2), and of dimension D = 3, namely 5* 3 , HP 3 , 
1/3,2 ~ 50(3), 1^3,3 ~ 0(3). We also made simulations for the D = 7 model V^. 

There exists an isomorphism between the manifolds 1RP 3 ~ V3 2 ~ 5*0(3). The Hamil- 
tonian are however different, except for the HP 3 model which can be mapped into the 
50(3) chiral model M. There are also the local isomorphisms between the manifolds 
5 1 ~ 0(2) and 5 3 ~ HP 3 ~ K 3 , 2 ~ 50(3) ~ 0(3). Concerning the global properties, 
the homotopy groups are n (O(2)) = n (O(3)) = Z 2 (boundary defect or domain walls), 
ni(]RP 3 ) = ILi(V 3 #) = Z 2 and Hi(5 x ) = Z (point defects). For the manifolds 5 3 and V 5 , 2 , 
the first homotopy groups n , Hi and n 2 are trivial (no topological defects of any kind, 
including instanton). 

Let us now turn to the RG (3 function for these models. For the S n ~ l models, we have 
the 5 loops results of Wegner J7J. As the manifolds ]RP n_1 are locally identical to 5 n_1 , 
the models have an identical (3 function (at all order), with the following relation between 
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the temperatures T(IRP n_1 ) = 2T(S' n ~ 1 ). For the 0(3) model, as the low temperature 
continuum limit neglects the two sheets of 0(n), it is identical to the SO (3) chiral model, 
and hence, via the previously mentioned mapping, to the HP 3 model. Therefore, the (3 
function is the S 3 (3 function (again, at all order), with the relation T(0(3)) = 8T(S 3 ). 
For the V n ^ models, the (3 functions have been computed by various authors |||8]|| and 
|T0| where the computations are more detailed. In the continuum limit, the action can be 
parametrized as 



f3H= — / d 2 x 

H 2T 1 + a 1 



(V<^) 2 + (V</? 2 ) 2 + {a- l)(<^iV^ 2 



(8) 



with the two coupling constants T and a. The bare models correspond to a = 1. At one 
loop order, the f3 functions in d = 2 + e are 

/l + a\ 2 T 2 

* = £ r-(n-2)(^) - (9) 

"""IK^) 2 *'"" 1 ^"*"" 3 ' 1 - (10) 

For d — 2, a trajectory of the flow possesses the invariant K = T\{n — l)a— (n— 3)\( n ~ 2 ^( n ~ 1 \ 
Then, we can integrate by the method of characteristics the flow equation and the differential 
equation for £. With the value a = 1 corresponding to the bare V n p models, we obtain 



Pi 



f3 a d a ln£ 



T 2 n — 2 

= — — r (11) 

27T -1 I n^l v 1 I -1 x 

a ~ 1 ~f n-2 ^k=0 ( n _i) fc+n _ 2 \ n -2 



By comparing eq. (|TT|) with eq.(|S|), we see that the term 1 — /? Q <9 Q ln£ changes only the 
proportionality constant in f3r- For n = 3, we have computed the above function at 2 loops 
order. In the domain of temperature we are interested in, the two loops correction is of 
order 2 ~ 3% and therefore can be neglected. For the S* 3 model, in the range of correlation 
length accessible to MC simulations, the 5 loops result gives a correction of order 10% over 
the one loop result. Similarly, for the MP 3 and 0(3) models, the 5 loops (3 function gives a 
correction of order 2%. Clearly, the lowest order f3 function is already a good estimate for 
all models. 

The Monte Carlo program has already been described in [[J. Briefly, we use the 
Ferrenberg-Swendsen multi-histograms method [II| and a cluster update algorithm [F2 



eventually biased according to |I3" |. The correlation length is computed based on the small 
k behavior of the Fourier transformed two points Green function G(k). From the histograms, 
we get an estimate for £(T) from which we compute l/c^rln£. Let us emphasize again that 
we need very good statistics in order to obtain a satisfactory estimate for this last quantity. 
The total computation time for the present simulations is within the range of several years on 
a typical RISC workstation. The largest size for which we can attain good enough statistics 
is typically L = 128. For the renormalization of the field rj(T), we can obtain a rough esti- 
mate, but our statistics are still insufficient to iron out the wrinkles on the curves. Finally, 
even with a cluster update, the autocorrelation time grows badly around the transition. 

Due to a lack of space, we can present the graphical results for only three representative 
models (by "representative" we mean that the f3 functions look similar). Let us emphasize 



5 



that there are no adjustable parameters in these results, neither in the MC measurements 
nor in the theoretical RG (3 functions. Figure 1 presents the results for the model, and 
is also representative of the S 3 model. The low temperature asymptotic regime is entered 
for correlation length of order £ ~ 10, and the measured and computed (3 functions match 
nicely. To the best of our knowledge, this is the first check of the 2 + e RG idea for a model 
different from the Heisenberg model. In the and S 3 models, there are no topological 
defects. 

Figure 2 presents the results for the HP 3 model and is also representative of the an d 
S 1 model. There we can see the clear discrepancy between the RG prediction and the MC 
measurement. The measured (3 function goes even below the RG prediction, a disturbing fact 
(although it does not violate any known rigorous bound). More disturbing is the fact that for 
T(1RP 3 ) ~ 0.3, the equivalent S 3 temperature is T(S 3 ) = T(IRP 3 )/2 ~ 0.15, a temperature 
way inside the asymptotic regime with a correlation length of order £(S 3 ,T = 0.15) ~ 10 8 . 
This series of models HP 3 , and S l contains point defects as IIi ^ 0. In particular, for 
the V 3:2 model used to describe frustrated antiferromagnets, this should be a clear warning 
(not to speak about the e = 1 extrapolation ||). 

Figure 3 presents the results for the 0(3) model, and is also representative of the 0(2) 
model. The discrepancy between the RG and MC (3 function appears even more obviously. 
Again, through the equivalence with S 3 in the continuous limit, the temperature T(0(3) ~ 1 
ends up way inside the asymptotic region for the S 3 model, also with a correlation length 
of order £ ~ 10 8 . The 0(2) and 0(3) models have boundary defects as n (O(n)) = Z 2 . 

From our grouping, the key factor appears to be the various topological defects. In our 
MC simulations, we also measured the density of the different defects and the transition is 
always clearly related to them. Obviously, there is a problem with the d = 2 RG (3 function 
for models with topological defects. In contradiction with the RG spin waves assumption, 
the global properties of the manifold do matter and a topological term seems to be a relevant 
direction for the RG flow. In fact, in view of the very different MC results for S 3 , MP 3 and 
0(3), it is even surprising that, for the S 3 model, the asymptotic prediction does not fail 
for £ ~ 10 8 but works down to correlation lengths of order £ ~ 10. 

Let us now speculate about the possible scenario in two dimensions. From the above 
figures, it is very tempting to conclude to a phase transition at a finite temperature for all 
models with topological defects. For models with boundary defects Ho ^ 0, similarly to 
the Ising model, we have a linear (3 function and a power law diverging correlation length. 
A finite size scaling for the correlation length and specific heat gives v = 0.78 for 0(2) 
and v = 0.71 for 0(3). Another possibility for this family of models would be to have 
two transitions, the one at a higher temperature corresponding to the symmetry breaking 
0(n) — > SO(n). As the broken symmetry is discrete, it can be accompanied with long range 
order. For models with point defects LTi ^ 0, similarly to the X-Y model, the (3 functions 
are compatible with a Kosterlitz-Thouless behavior (3 ~ \T — T c \ 3 / 2 . Let us emphasize that, 
for model with topological defects, phase transitions at a finite temperature are completely 
compatible with a finite size scaling analysis of the data. In contradistinction, the above RG 
predictions are consistent with the simulation for models with trivial topologies, but clearly 
ruled out for model with topological defects. Of course, the MC simulations and finite size 
scaling are always powerless against arguments of the kind: for a correlation length of order 
£ ~ 100 (or £ ~ 10 12 ) there is a crossover to the predicted RG asymptotic behavior. What 
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is missing in order to draw a firmer conclusion is a theoretical basis. We would like for 
example to have a modification of the Villain mapping for models with curvature. It is easy 
to build an energy-entropy argument for a pair of point defects, but it seems difficult to 
turn this simple argument into a good approximation. For a model believed to be in the 
same universality class as the model, from an energy-entropy argument and from MC 
simulation, Kawamura and Miyashita |14| already proposed a phase transition at a finite 
temperature. Beside possibly dangerous analogies with the X-Y model, only the large n- 
limit for IRP n concludes to a (first order) transition at a finite temperature . Caracciolo 



et all [IB], studying a mixed S 3 — IRP 3 model, concluded also that the usual RG computation 
does not describe the HP 3 model and proposed some interesting alternatives, although still 
favoring a transition at T c = for the pure HP 3 model. Another possible scenario would 
be that the naive power counting used to compute the relevance of various term is valid 
only very close to T = 0. At higher temperature, the present results indicate that the least 
irrelevant terms are the one related to topological defects. 

Simply stated, in the range of temperatures accessible to MC simulations and by using 
finite size scaling, our data for all models with topological defects are very well described by 
a phase transition at a finite temperature. 

I would like to thanks B. Halperin, D. Nelson and U. Tauber for various discussion. This 
work is supported by the Fonds National Suisse de la Recherche Scientifique, by the NSF 
grant DMR-91-15491 and by the office of Naval Research grant N00014-93- 1-0190. 
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FIGURES CAPTIONS 

Figure 1: The theoretical (dotted curve) and measured (3 function for the V 5>2 model. 

Figure 2: As figure 1, but for the IRP 3 model. At this scale, the one and five loops 
theoretical results are undistinguishable. 

Figure 3: As figure 2, but for the 0(3) model. 
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